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Abstract: Liquid xenon detectors such as XENON10 and XENON100 obtain a significant fraction 
of their sensitivity to light (< 10 GeV) particle dark matter by looking for nuclear recoils of only 
a few keV, just above the detector threshold. Yet in this energy regime a correct treatment of 
the detector threshold and resolution remains unclear. The energy dependence of the scintillation 
yield of liquid xenon for nuclear recoils also bears heavily on detector sensitivity, yet numerous 
measurements have not succeeded in obtaining concordant results. In this article we show that the 
ratio of detected ionization to scintillation can be leveraged to constrain the scintillation yield. We 
also present a rigorous treatment of liquid xenon detector threshold and energy resolution. Notably, 
the effective energy resolution differs significantly from a simple Poisson distribution. We conclude 
with a calculation of dark matter exclusion limits, and show that existing data from liquid xenon 
detectors strongly constrain recent interpretations of light dark matter. 
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1. Introduction 

The recent results of an 11 live day direct search for dark matter from the XENON100 collaboration 
PJ have brought liquid xenon detectors into sharp focus. If a dark matter particle scatters in 
the 40 kg active target, the experiment hopes to detect the scintillation photons and ionisation 
electrons that are produced by the recoiling xenon nucleus. This class of technology, which notably 
also includes the next-generation LUX experiment [g, [| , has shown excellent scalability: the 40 kg 
fiducial (out of 161 kg total) 11 live days exposure from XENON100 offers a similar sensitivty to 
the approximately 128 live days exposure from the CDMS collaboration JO, due almost entirely 
to the increased target mass. For reference, the state of the art in liquid xenon detectors in early 
2007 was a 7.2 kg fiducial target (out of 31 kg total) ||. However, there is significant disagreement 
regarding the sensitivity of liquid xenon detectors to light (< 10 GeV) dark matter particles |J 0, |§| . 
The argument centers on the non- linear scintillation yield of liquid xenon for nuclear recoils, C e ff, 
and in particular on extrapolations of £ e ff below the lowest energy at which measurements are 



available. Since liquid xenon dark matter experiments use the primary scintillation signal to set 
the energy scale for particle interactions, uncertainty in C e ff leads to uncertainty in the deposited 
recoil energy. Near the detector threshold of a few keV, the uncertainty is more properly a question 
of whether or not a low-energy nuclear recoil offers any detectable signal. For light dark matter, 
this can have order-of-magnitude implications for which values of mass m x and cross section a x of 
a hypothetical dark matter particle are consistent with the results of the experiment. 

The energy dependence of C e ff nas been measured by numerous groups 0, [HI 11, 11, 13, 14 



and substantial disagreement persists — particularly at nuclear recoil energies E„ r < 10 keV. The 
measurement is very challenging in this regime: a 5 keV nuclear recoil is only expected to yield 
1.6 (2.4) detectable scintillation photons in XENON100 (XENON10), if £ e // = 0.09 at that energy 
|l|, P-5| . It would seem that an alternate measurement technique is needed, or failing that, an 
additional means to constrain the measurement. Fortunately, such a means exists, although it 
has not previously been exploited for this purpose. A particle interaction in liquid xenon creates 
a number N e of ionisation electrons in addition to iV 7 scintillation photons. The ratio N e /N 1 
is already used to discriminate between nuclear recoils (from neutrons, and expected from dark 
matter) and electron recoils (from background gamma and beta radiation). 

In Sec. we show how N e /N-y can be used to constrain C e ff, by requiring consistency with 
measurements of N e /N 7 as a function of energy. We show reasonable lower, central and upper 
values of £ e // in Fig |[ It turns out that some extrapolation of £ e // is necessary and justified, even 
below the lowest measured experimental values; this is also shown in Fig. ra. A correct calculation 
of dark matter sensitivity requires a correct treatment of the low-energy detector efficiency, which is 
discussed in Sec. |l]and Sec. [h| for XENON10, and in Sec. | for XENON100. Although Poisson 



fluctuations are the dominant component of the energy resolution, the effective energy resolution 



departs significantly from the Poisson prediction. This is discussed in Sec. 4.3, and is due in large 
part to the fact that experiments only look for dark matter signals in the lower half of the nuclear 
recoil band. In Sec. @ we apply all these results to a calculation of dark matter exclusion limits. 

2. Preliminaries 

2.1 Signal production and measurement in liquid xenon 

This discussion is applicable to the XENON10, XENON100, Zcplin II [§, Zcplin IIlQ and LUX 
experiments. Numerical examples will use properties of the XENON10 detector |fl6| . A particle 
interaction in liquid xenon creates both excited Xe* and ionised Xe + xenon atoms [|17| , which react 
with the surrounding xenon atoms to form excimers. The excimers relax on a scale of 10 -8 s with 
the release of scintillation photons. This prompt scintillation light is detected by photomultiplier 
tubes and is referred to as the SI signal. The total efficiency for a single photon to produce a 
single photoelectron in a photomultiplier tube is ~ 0.10 in XENON10, and varies by detector. It 
depends primarily on the product of the geometric light collection efficiency of the detector, and 
the quantum efficiency of the photomultipliers llqj . 

An external electric field Ed across the liquid xenon target causes a large fraction of ionisation 
electrons to be drifted away from an interaction site. The electrons are extracted from the liquid 
xenon and accelerated through a few mm of xenon gas by a stronger (by about x 10) electric 
field E e , creating a secondary scintillation signal. This scintillation light is detected by the same 
photomultiplier tubes and is proportional to the number of ionisation electrons. It is generally 
referred to as the S2 signal. XENON10 measured approximately 25 S2 photoelectrons per extracted 
electron @. 

Dark matter search experiments have tended to report the response of the detector to particle 
interactions in terms of logio(S2/Sl) as a function of SI fy, g, [n], [15| Q. When such a plot is 



produced as a result of a dedicated neutron (gamma) calibration, the result is generally referred to 
as the nuclear (electron) recoil band. Typical examples can be found in Jj], Q. For nuclear recoils, 
the measured quantity SI (in units of photoelectrons) is related to the nuclear recoil energy E nr via 

' S ' S ' (2.1) 



For XENON10, the light yield L y = 3.0 photoelectrons/keVee 1 for 122 keV photons, and the 
scintillation quenching of electron and nuclear recoils due to the electric field E^ = 0.73 kV/cm are 
S e — 0.54 |14J and S n — 0.95 ||. These quantities depend on the applied electric field, which is 
similar for XENON100 (E d = 0.53 kV/cm) and higher for ZEPLIN III (E d = 3.9 kV/cm). This 
construction of C e ff in terms of SI and the 'standard candle' L y is necessary because the probability 
of detecting a single scintillation photon is significantly less than one. Physically, L y /S e represents 
the total detection efficiency for primary scintillation photons. 

In contrast, the measured quantity S2 (also in units of photoelectrons) is related to the nuclear 
recoil energy E„ r via 

E„ r = — , (2.2) 

where Q y is the number of detected electrons per keV. It is possible to express Q y in absolute terms 
because single electrons are measured with high efficiency and relative ease. Both of these quantities 
are reported in [|9| , in which the true nuclear recoil energy is known from the initial neutron energy 
and the scattering angle. Note that in contrast to C e ff, Qy is not referred to zero electric field. 
This may lead to some residual dependence on the electric field. 

2.2 Monte carlo simulation of signal production and measurement 

Because Q y and C e ff govern signal production from nuclear recoils in liquid xenon, they must also 
dictate the shape of the nuclear recoil band. This observation, which has not previously been used 
to constrain measured values of Q y and £ e ff, is the premise of Sec. |3J. Given Q y and £ e ff as a 
function of nuclear recoil energy, it is relatively easy to generate a nuclear recoil data set via monte 
carlo simulation. The results of this simulation will be compared with the measured XENON10 
nuclear recoil band |Tq], for various combinations of Q y and £ e ff- In the interest of clarity, the 
following description of the simulation will therefore use numerical examples from the XENON10 
detector @. 

For each simulated event with nuclear recoil energy E nr , which was modeled in 0.25 keV steps, 
the simulation output is a value of SI and S2 as would be seen by the detector. This was accom- 
plished as follows: given E nr , the expected number of detected electrons N e was calculated from Q y 
and subjected to Poisson fluctuations. The number of S2 photoelectrons was then obtained from 
N52 = N e • 25 photoelectrons/electron, taking care to account for the measured width (a/ /J, = 0.30) 
of the single electron distribution ]lq| . The expected number of SI photoelectrons Ngi was cal- 
culated from C e ff, modeled as a discrete random hit pattern on the photomultipliers, and also 
subjected to Poisson fluctuations 2 . Both SI and S2 were subjected to Gaussian fluctuations in the 
size of a single photoelectron measured by the photomultipliers, using the average a j ' jjl = 0.58 |16| . 
Finally, the SI signal was required to satisfy ann>2 coincidence requirement in the photomulti- 
pliers [n9| , within a 300 ns time window |12| . The primary scintillation signal was assumed to have 
a pulse shape as shown in Fig. [LOL The scintillation decay time is discussed in the Appendix. 

lr The suffix 'ee' refers to electron recoil equivalent energy. The suffix 'r' is usually appended to indicate nuclear 
recoil equivalent energy, to emphasize the quenching of the electronic signal from nuclear recoils. In this work we 
always mean nuclear recoil energy unless stated otherwise, and will reserve the 'keVr' unit to explicitly indicate 



measured nuclear recoil energy. This is discussed in detail in Sec. 4.1. 

2 Strictly speaking, the probability for a scintillation photon to produce a photoelectron is governed by binomial 
statistics. In practice, Poisson statistics are a very good approximation. 



It is well-known that the electron and photon signals in liquid xenon are anti-correlated, for 
electromagnetic scattering [g0| glj. Anti-correlation specifically refers to the situation in which an 
ionisation electron recombines with a XeJ , leading to the release of a scintillation photon. This 
is described in detail by M Q4| and others. Recombination fluctuations are evident in gamma 
lines |12L EU|, and contribute significantly to the observed width of the electron recoil band pL 19 



Surprisingly, there is scant evidence of significant recombination fluctuations for nuclear recoils. 
Recent modeling p2] of signals from nuclear recoils in liquid xenon explains this as originating 
from an initial production ratio Xe*/Xe + ~ 0.9 for nuclear recoils, compared with Xe*/Xe + = 
0.06 for electron recoils. Physically this means that for nuclear recoils, much less of the recoil 
energy produces ionization; this in turn significantly limits any potential effect of recombination 
fluctuations. It may also explain why the nuclear recoil band falls below the electron recoil band. 
The initial production ratio of Xe*/Xe + is weakly dependent on Ed, considering the available data 
on scintillation quenching 0, H . There may also be a small dependence on nuclear recoil energy, 
due to track structure of the recoiling nucleus. We assume that these details are already encoded 
in the observed nuclear recoil band iV e /AT 7 , reported as logio(S2/Sl). 

Effects arising from recombination fluctuations were deemed to be sub-dominant for nuclear 
recoils and were not modeled in this simulation. If this assumption were incorrect, the simulated 
width of the nuclear recoil band would be too small. As we shall see in Sec. |3j, the simulation is 
able to make a very reasonable reconstruction of the observed nuclear recoil band width. A note 
on boundary conditions is in order. In all cases, we set Q y = and C e ff = at E nr < 0.5 keV. 
However, in all cases, the SI n > 2 coincidence requirement rejects events with such low energy. 
This will be discussed in more detail in Sec. 0. 



3. Nuclear recoil band constraints on Q y and C e ff 

We will consider three case studies to illustrate the utility and limitations of using the nuclear recoil 
band to constrain Q y and C e ff- This is simply a statement that the energy dependence of N e /N 1 , 
and therefore also Q y /C e ff, are encoded in the reported nuclear recoil band logio(S2/Sl). The 
cases are shown in Fig. p] as solid curves, each of which was modeled as a cubic spline. They are 
truncated at the value of E nr below which no signal would have been recorded. No extrapolation 
of Q y or C e ff was made (nor would it be relevant), beyond the extent of the curves shown in Fig. 
|l|. This cutoff is enforced primarily by the SI n > 2 coincidence requirement, in conjunction with 
the 12 electron S2 threshold. 

3.1 Case 1 

As a starting point, we take a spline through the central values of Q y and C e ff from S and check if 
our simulation reproduces the measured XENON10 nuclear recoil band (which we take from p5|) . 
These curves are labeled Case 1 in Fig. |l|. A notable difference between the Manzur et al., M 
experiment and |L5| is the different electric field Ed across the active target. While the difference 
between Ed — I kV/cm and Ed = 0.73 kV/cm appears to be negligible for scintillation quenching 



of nuclear recoils [h4l, it may have a small effect on Q y . We will discuss this further in Sec. 3.2. 

Figure || shows the XENON10 nuclear recoil band /i and ±3<r contours (red stair steps), along 
with the 12 electron S2 threshold (solid red curve) described in p5||. The simulated data points 
are shown in light gray, and data points corresponding to E nr = 5 keV and E nr = 30 keV are 
highlighted in dark blue to elucidate the mapping between the nuclear recoil band and true nuclear 
recoil energy. The XENON10 collaboration report that their nuclear recoil band was calculated in 
discreet SI bins, by fitting a Gaussian to the logio(S2/Sl) distribution |L5|. This distribution is 
empirically highly Gaussian fig], with some non-Gaussian tails evident on the low side. 
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Figure 1: Q y and £ e // values used as input to the simulation (solid curves), for each case study. Also 
shown are data from |] (blue •, E d = 1.0 kV/cm), @ (0), jjj) (V, E d = 0.73 kV/cm) and §|] (°. 
Ed = 0.1 kV/cm and ■ Ed = 2.0 kV/cm). The top panel also shows the trend of Q y corresponding 
a measurement of 30 electrons, regardless of nuclear recoil energy (red dash). Q y may depend on E e in 
addition to the weak dependence on Ed; the former is often not explicitly reported. 
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Figure 2: The XENON10 nuclear recoil band parameterization centroid and ±3a (red stair steps) |l5[| are 
compared with the simulated nuclear recoil band data (light grey points). Simulated data corresponding to 
E nr = 5 keV and E nr = 30 keV are highlighted in dark blue to indicate the mapping between the nuclear 
recoil band and true nuclear recoil energy. Resulting centroid and ±3<r contours from Gaussian fits to the 
simulated data are indicated in the figure legend. The S2 software threshold that XENON10 implement is 
duplicated in the simulation. 



We employed exactly the same procedure and bin widths in calculating the nuclear recoil band 
obtained from our simulation, and display the /i (as o) and ±3cr (as A and V) values, with horizontal 
bars indicating the bin width. Simulated events which appear below the 12 electron threshold were 
discounted prior to making a Gaussian fit. The statistical error in the simulated fj, and ±3ct values, 
if not visible, is smaller than the size of the data points. Bearing in mind that the y axis in Fig. |2| 



is logarithmic in S2/S1, the systematic disagreement between our simulated fi and the XENON10 
measured fi is severe. Evidently, self-consistency requires a decrease in Q y , or an increase in C e ff. 
Although it is not discussed in detail here, we also performed this same analysis using a spline 
through the central Q y and C e ff values reported by XENON10 Q. These values are indicated 
by V in Fig. El The simulated fi were systematically lower than the measured XENON10 fi, 
indicating that either an increase in the measured Q y , or a decrease in C e ff would be necessary 
for HH to obtain self-consistency. An example of the requisite decrease in £ e // that would yield 
self-consistency is provided by the dotted curve in Fig. @. 

3.2 Case 2 

The Q y indicated by the curve labeled Case 2 in Fig. fil was constructed to obtain consistency 
between the central values of C e ff reported in M, the XENON10 nuclear recoil band and the 
simulated nuclear recoil band. The agreement is quite good, as shown in Fig. [|. Note that this 
Q y curve is better than la consistent with the measurement of ||, except below about 8 keV. The 
fact that the XENON10 data in this case require a systematically lower Q y at all energies may 
simply be a reflection of the differences in Ed, as discussed in Sec. 3.1. The flat extrapolation of 
C e ff — 0.07 between 1.2 keV and 4 keV was necessary to align the fj, values of the nuclear recoil 
band, in the 2 — 5 SI photoelectron bin. The falling extrapolation of £ e ff as illustrated by Case 1 is 
perhaps more aligned with theoretical expectations [24]. However, agreement with the XENON10 
nuclear recoil band in that case would have required Q y to be 2<r below the measured values B at 
4 keV and 5 keV. In general one should expect a measurement of Q y to be significantly more robust 
than a measurement of £ e //. For example, a typical 6 keV nuclear recoil in the XENON10 detector 
has the expectation (from Fig. ||) of 3 SI photoelectrons and about 750 S2 photoelectrons (30 
ionisation electrons). In contrast to the SI signal, the S2 signal is (in general) far above threshold 
and far from the regime in which Poisson fluctuations dominate. These considerations imply that a 
Q y curve lower than indicated by Case 2 would be difficult to explain. As it stands, there remains 
some tension with the measured value at 4 keV ||. This would argue for higher values of C e ff 
below 8 keV. 
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Figure 3: All symbols and curves as described in Fig. 



A possible resolution lies in the S2 trigger threshold. The detector used in || appears to have 
an S2 gain of about 8 photoelectrons per ionisation electron (about x3 lower than XENON10), 



based on the mean value of ~ 300 S2 photoelectrons measured for E„ r ~ 6 keV (see [[|, Fig. 8c). 
Suppose that the S2 trigger efficiency begins to drop below about 180 S2 photoelectrons (as it does 
in |16J ) . This would cut off the low side of the measured S2 distribution, and thereby artificially 
inflate the measured central S2 value at low recoil energies. Figure [J shows the Q y curve that would 
result if fl measured a central value of 240 S2 photoelectrons (30 electrons) at all recoil energies 
(red dash curve). Above E„ r ~ 8 keV this is clearly unrealistic, but for the lowest nuclear recoil 
energies it could indicate a measurement trend imposed by the S2 trigger efficiency of the detector. 
This could ease the tension between the measured value of Q y at 4 keV, and the Q y curve required 
for simultaneous consistency with the central £ e ff values of || and the XENON10 nuclear recoil 
band. 

3.3 Case 3 

In order for the Zeplin III measurement of C e ff to achieve consistency with the XENON10 nuclear 
recoil band, an abrupt drop in Q y is required. The level of consistency is not shown, but appears 
similar to that shown in Fig. (without the dark blue band corresponding to 5 keV events). The 
requisite crash in Q y , labeled Case 5 in Fig. [j], cannot be reconciled with measurements J9J. Given 
the robustness of the S2 signal, we are forced to consider other explanations. One possibility is that 
the electron extraction efficiency in the Zeplin III detector is significantly lower than in XENON10. 
This effect is expected, although there is some uncertainty in its magnitude, because the electric 
field across the liquid-gas interface is about 25% lower in Zeplin III than in XENON10. However, 
such an effect would be expected to result in a uniform systematic decrease of Q y as a function of 
recoil energy. Other possibilities are that the Zeplin III SI detection efficiency may be lower than 
stated in Fig. 14 of 0. 

3.4 Summary: most reasonable estimate of Q y and C e ff 

At this point we pause to make a few observations. First, the Q y and C e ff curves labeled Case 2 in 
Fig. fy are reasonably consistent with recent measurements M, [Hi E3l, and the simulation indicates 
consistency with the XENON10 nuclear recoil band fj, and a. The reproduction of the observed 
XENON10 band width a confirms our assumption that recombination fluctuations are subdominant. 
The Case 2 Q y curve is as low as can be comfortably explained by a possible systematic effect 
between the slightly different electric fields Ed, given the near absence of scintillation quenching for 
nuclear recoils fbl. It is important to recognize the most glaring weakness of this analysis: the lack 
of a complete set of Q y , £ e ff and nuclear recoil band measurements under the same conditions. 



This inevitably leads to systematic uncertainty. As mentioned in Sec. 3.1, the magnitude of this 
uncertainty does appear to be small. We therefore propose that the central values of ||, with 
the extrapolation indicated by Case 2, represent the most reasonable estimate of C e ff given the 
available data. 

3.5 A proposed calibration strategy 

Because of the tendency of liquid xenon detector experiments to report their results with an Sl- 
based energy scale, an enormous amount of pressure has been placed on measuring C e ff- This has 
shifted the experimental focus away from what should be a far easier measurement, that of Q y . The 
results of Sec. suggest a clear course of action for future tagged neutron scattering measurements 
of C e ff, as there will undoubtably be. In addition to the methods and techniques already employed 
p, 10 an additional first step should be the acquisition of nuclear recoil data with a broad-spectrum 



neutron source such as americium-beryllium, or californium. This will provide a nuclear recoil band 
and an essential cross-check. The second step should be to make a robust measurement of Q y using 
a tagged neutron source of known energy. In order to minimize the measured uncertainty on Q y , 



the gain of the S2 signal should be high, > 25 photoelectrons per ionisation electron. Ideally the S2 
trigger threshold would be set at the level of a single ionisation electron, as XENON10 have proven 
possible |Ea] . If that were to lead to an excessive trigger rate, a careful assessment of the trigger 
threshold with a voltage step input [161 could be used instead. This, combined with the analysis 
detailed in Sec. 0, should allow a significant reduction in systematic and statistical uncertainty in 
the energy dependence of C e ff. 

4. Low-energy nuclear recoils 

The discussion in Sec. [| assumes Q y and £ e // indicated by Case 2 in Fig. [EL and will continue to 
assume instrumental parameters corresponding to the XENON10 detector. 

4.1 Detector efficiency 

Event detection efficiency r\ as E nr approaches zero is dominated the SI detection efficiency. An 
S2 software threshold at the level of 12 electrons (XENON10) or 15 electrons (XENON100) plays 
a very minor role. Our simulation returns a direct measurement of 77, which is shown in Fig. ( V) 
as a function of kcV. 
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Figure 4: The detection efficiency »?(Sl) (dotted curve), ry(keVr) (dashed curve) obtained from Eq. 2A, 
ry(keV) (solid curve) obtained from Eq. j3| and r\ as obtained directly from the simulation (V). Also shown 
are the acceptance of the '50% acceptance box' (a) and the combined low-energy acceptance for nuclear 
recoils (0). The x axis is in keV unless indicated otherwise in the legend. Statistical uncertainty, where not 
visible, is smaller than the data point. 



A more accessible quantity is frequently reported, the behavior of 7?(S1) g, [l(], |l2[ Q. It 
can be obtained by a simple photoelectron statistics monte carlo simulation, and is a function 
of the number of detected SI photoelectrons. Comparison with low-energy Compton scattering 
background (e.g. Fig. 4.8 of |l"|) can be used to confirm its validity. In order to correctly apply 
the detector efficiency 77 to the application of dark matter exclusion limits, it is necessary to obtain 
either ?7(keV) or 77(keVr). The distinction is important here. As previously mentioned, keVr is 
already widely used to mean 'keV nuclear recoil equivalent energy,' to emphasize the quenching of 
the electronic signal from nuclear recoils. The keVr unit carries the tacit assumption that it is a 
measured quantity. It must therefore also carry the imprint of the detector resolution. In contrast, 



the keV unit is most apt for a calculated or theoretical quantity, with perfect resolution. There are 
two correct paths forward: 

1. convert //(SI) to ?7(keVr) via C e ff- The result should only be applied to a theoretical prediction 



of the differential scattering rate, such as Eq. 6.1, after first correcting for the detector 
resolution. 



2. convert ry(Sl) to 7?(keV) via the method described in Sec. 1.1.2 . The result can then be 
applied directly to a theoretical prediction of the differential scattering rate, since account 
has been taken of detector resolution. 

4.1.1 How to obtain 77(keVr) 

In Fig. [| we show r?(Sl) for Case 2 (dotted curve). It agrees quite well with 7?(S1) as reported 3 by 



the XENON10 collaboration jl2[. As an example, consider that 77(2. 1)= 0.8. From Eq. 2.1 and the 
C e ff curve represented by Case 2, we obtain from the conversion 2.1 SI photoelectrons = 4.7 keVr. 
Thus, 77(4.7 keVr) = 0.8. Applying this procedure to the entire ?7(S1) curve gives 7i(keVr), as shown 
in Fig. (dashed curve). While both 77(S1) and ?y(keVr) are correct representations of the detector 
efficiency, neither can be applied directly to a theoretical prediction of the differential dark matter 



scattering rate (such as Eq. 6.1). The difference between these two curves and ry(keV) lies in the 
detector resolution. 

4.1.2 How to obtain 77(keV) 

To be explicit, £ e ff is a measure of the expected number of detected SI photoelectrons, given 
a nuclear recoil energy E nr (in keV). In calculating the nuclear recoil energy E„ r , given some 
number of detected SI photoelectrons, the correct response is to specify a value of E nr and a 
probability that it resulted in the observed SI measurement. As an example consider a measurement 
of 2 SI photoelectrons. Such an event is the expected result for an initial E„ r = 4.6 keV, but 
E„ r = 8 keV still has a 10% chance of resulting in SI = 2 photoelectrons. The probabilities are 
governed by £ e // and Poisson fluctuations. To correctly calculate 77(4.6 keV), we must account for 
all possibilities. We do this by piece-wise multiplication of 7y(Sl) with the Poisson probabilities that 
a measurement returned [0123 ...] photoelectrons, given the expectation k = 2 photoelectrons. In 
other words, define 

r]f(k) = fpoi SS on([0 12 3...];k) (4.1) 

for integer values of k detected photoelectrons. For k = 2, 

fpotssondO 1 2 3 ...]; 2) = [0.14 0.27 0.27 0.18 ...]. (4.2) 

The total detection efficiency is given by 

Tikevity^^Vfi^^-ifsiii) (4.3) 

i 

in which the nuclear recoil energy E„ r is obtained for each k via C e ff- For k = 2, we obtain 



77(4.7 keV)= 0.5. The predicted 77 obtained from Eq. 4.3 is shown in Fig. (solid curve), and 



agrees well with 77 obtained directly from the simulation (v). The domain of each 77 appearing in 



Eq. 4.3 is explicitly subscripted for clarity. Note that 7i(keV) is only defined for integer values of k 
detected photoelectrons. The solid curve in Fig. || truncates at E„ r ~ 2.8 keV, which corresponds 
to ft = 1 SI photoelectron. 
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Note that in Fig. 5 of [113] , r)(Sl) is referred to as 'SI peak- finding,' and is shown with a factor x3 in the x axis. 



4.2 The 50% acceptance box 

Experiments |y, |l5[ [19| have tended to use the region of the nuclear recoil band between \x and 
fi — 3cr as the dark matter "50% acceptance box." Considering the tendency of experiments to 
report results as shown in Fig. pi one should expect 50% acceptance as a function of SI. In 
Fig. f| we show the acceptance of the 50% acceptance box for nuclear recoils (A) as a function 
of keV. For E„ r > 20 keV, the 50% acceptance box lives up to its name; as E nr approaches zero, 
the acceptance of the region between /z and /i — 3a climbs by up to 30%. This should not be 
surprising: at small nuclear recoil energies, the only events which experiments are able to record 
have benefitted from upward Poisson fluctuations in the number of detected SI photoelectrons. 
Assuming normally distributed S2 values 4 , such events are forced to lower values of logio(S2/Sl), 
and so in a typical nuclear recoil band plot (as a function of SI), appear preferentially below the 
measured band centroid. The acceptance (a) shows more scatter than might be expected from the 
stated uncertainty; this is simply an effect of the discreet SI bins employed in the calculation. The 
combined low-energy acceptance of the 50% acceptance box and r] is also shown in Fig. || (0). 

Although the nuclear recoil band as a function of detected SI photoelectrons is empirically 
highly Gaussian, there is a small non-Gaussian component at low recoil energies. This non-Gaussian 
component is a natural result of Poisson fluctuations, although other factors that our simulation 
does not consider may also contribute. 

4.3 Detector resolution 

As mentioned, experiments have tended to report event energies in terms of detected SI photo- 
electrons p], 01| fL9[ , One could reasonably expect the detector resolution for nuclear recoils, as a 
function of E nr , to be governed by Poisson statistics for the number of detected SI photoelectrons. 
Naturally, reality is more complicated. 
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Figure 5: (top) SI distribution for E nr =5 keV, and Poisson distribution (blue dash) for the expected 
number (k = 2.3) of SI photoelectrons at this recoil energy, (bottom) The same SI distribution, and 
Poisson distribution for k — 2.3, both transformed to keVr via £ e // from Fig. bl A Gaussian fit (red dash 
dot) is also shown. 



The simulation allows us to directly measure the detector resolution, considering the inputs 



listed in Sec. 2.2. The distribution of SI photoelectrons for events falling in the 50% acceptance 

4 This is reasonable; from Fig. pi even at E nr = 2 keV, one expects to detect ~ 12 electrons. 
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box is shown in Fig. [| for the case E„ r = 5 keV. Considering Case 2 shown in Fig. [y, XENON10 
should expect to measure 2.3 SI photoelectrons at E„ r = 5 keV. A Poisson distribution with k = 2.3 
is shown (top panel, blue dash), normalized to the total N = 5000 simulated events. The same 
distribution is shown (bottom panel, blue dash) transposed from SI photoelectrons to keVr using 
C e ff, and also fails to match the simulation result. A total of 1887 events comprise the histograms 
in Fig. pi This is reflected in the total acceptance (0) of 0.37 shown in Fig. ||. An additional 
1158 events with E„ r = 5 keV fell in the region between fj, and [i + 3cr, which is reflected in the 
acceptance (a) of 0.61 shown in Fig. ||. The remaining 1955 events either failed to satisfy the n > 2 
SI coincidence requirement, or the S2 > 300 photoelectron requirement. This is reflected in the 
acceptance (V) of 0.61 shown in Fig. [|. 

In terms of either SI photoelectrons or keVr, the observed distribution is at best quasi Poisson. 
The reason for this is two- fold: the left edge of the distribution is truncated by both the SI detection 
efficiency and the top of the 50% acceptance box. Both edges are broadened by photomultiplier 
fluctuations. Figure ra illustrates that at low recoil energies, a detector will retain sensitivity to 
upward fluctuations in measuring the primary scintillation signal. However, it also illustrates that 
one must be explicit about the underlying assumptions when "taking into account an SI resolution 
dominated by Poisson fluctuations" [EJ. For reference Fig. (bottom) also shows a Gaussian fit 
(red dash-dot), although the quality of the fit is rather poor. 
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Figure 6: All symbols and curves as described in Fig. B. The Poisson distribution (blue dash) has k — 21. 

A second example is shown in Fig. pi for E„ r = 30 keV. The left edge of the distribution is now 
truncated solely by the top of the 50% acceptance box. The right hand edge is visibly widened, 
relative to the Poisson expectation (blue dash), by photomultiplier fluctuations. In terms of keVr, a 
Gaussian fit (red dash-dot) provides a reasonable parameterization. This is true for E„ r > 25 keV. 

It seems unlikely that the detector resolution can be summarized with a simple statistical model, 
particularly at low recoil energies where it is most relevant. For light mass (< 10 GeV) particle 
dark matter, in which case detector sensitivity is dominated by Poisson fluctuations above the SI 
threshold, a correct model of the detector resolution is essential. As a practical matter, it seems 
preferable to use the discrete resolution distributions obtained directly from the simulation |26j. 



4.4 Cut acceptance e 

The cut acceptances e reported by e.g. the XENON10 collaboration ||, |j| or XENON100 collab- 
oration jy as a function of detected SI photoelectrons reflect the fraction of nuclear recoil events 
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passing all cuts, in each SI bin. A correct mapping between SI and keV requires considering the 
probabilities that a measured SI value corresponded to each value of keV nuclear recoil energy. In 
calculating dark matter exclusion limits, it is therefore more natural to scale the bin edges for each 
e to keVr via C e ff, they should then be applied to the theoretical prediction for the differential 
scattering rate (Eq. |6.1[) only after correcting for detector resolution. 



5. XENON100 

The results obtained in Sec. [| and Sec. |] also apply to the XENON100 experiment, with 
slight numerical modification. The major difference from a sensitivity point of view is the lower 
Ly/S e = 2.2/0.58 reported by the XENON100 collaboration Q. This is the primary reason that 
XENON100 has a reduced sensitivity to low mass dark matter, relative to XENON10. The su- 
perior sensitivity of XENON10 in the low mass range has already been noted |g7|], but attributed 
solely to the lower SI threshold (2 photoelectrons at XENON10, versus a stated 4 photoelectrons 
at XENON100). XENON100 reports @ a slightly lower S2 gain of about 20 photoelectrons per 
extracted electron. Practically, this appears to translate into a 15 electron (rather than 12 elec- 
tron) S2 software threshold. However, the S2 threshold has a minimal impact on rj, since 77 is still 
dominated by the SI n > 2 coincidence requirement and L y /S e . 




Figure 7: XENON100 detector acceptance, assuming Q y and £ e ff given by Case 2 from Fig. 
symbols and curves are as described in Fig. v&. 



All 



Finally, we note that the XENON100 collaboration have opted to require the SI n > 2 coin- 
cidence within a 20 ns time window m. There is some uncertainty in the scintillation decay time 
for nuclear recoils in liquid xenon, and in the Appendix we describe a reasonable parameterization. 
With this parameterization, we find that a 20 ns coincidence window (rather than 300 ns) reduces 
f] by 20% at 5 keV, and about 10% at 10 keV. In Fig. we show the detector efficiency obtained 
with the XENON100 detector parameters, given the physical assumptions of Case 2 from Fig. [j]. 
Excellent agreement with the XENON10 nuclear recoil band was obtained (but is not shown). This 
should not be surprising given the instrumental similarity between the two experiments. We note 
that the XENON100 collaboration have already confirmed the similarity of the nuclear recoil bands 
|Pq| obtained from the two experiments. The only notable difference we observe is a slight upturn 
in the simulated XENON100 \i values, below about 15 SI photoelectrons. This was also observed 
in the XENON100 data H. 
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6. Dark matter exclusion limits 

Having obtained significant insight on the low energy acceptance rj, the actual acceptance of the 
50% acceptance box and the detector resolution, we turn to calculation of dark matter exclusion 
limits. As in Sec. ra and Sec. m we use the Q y and £ e ff curves given by Case 2 as a central value. 
Noting that any reasonable C e ff curve must also lead to reasonable values for Q y , we choose the 
set of curves shown in Fig. || (dashed) as a lower bound. This is conservative considering that lower 
values of Q y would be very difficult to reconcile with measurements ||. We also consider a spline 
through the measured values of p| as an upper bound for £ e // , as shown in Fig. [| (dotted) . While 
higher values of C e ff could also be considered reasonable, we find it more interesting to consider 
a case which coincides with a published measurement. Each of these cases results in excellent 
agreement with the XENON10 nuclear recoil band, similar to that shown in Fig. ||. 




50 100 



Figure 8: Q y and £ e // values used to calculate central (solid curve, also shown in Fig. hi as Case 2) most 
conservative (dashed) and upper (dotted) 90% C.L. bounds. Each curve is extrapolated as necessary to 
reproduce the measured nuclear recoil band shown in Fig. pi The dotted curve is truncated at 1 keV. Other 
data as indicated in Fig. hi 



We calculate predicted differential event rates as a function of nuclear recoil energy on a xenon 
target from 

dR " — Px a n , 2 p2^ N [ V ° SC f( V ) 



= N T M^ 



n: 



-A'F z (E nr ) 



-dv, 



(6.1) 



where the number of target nuclei in the detector is Nx, the mass of the target nucleus is M^ and 
its atomic number is A. We assume the standard local dark matter density p x — 0.3 GeV cm~ 3 , 
with dark matter particle mass m x and cross section (per nucleon) a n . The reduced mass \x ne is 
for the nucleon — dark matter particle system. The nuclear form factor F(E nr ) accounts for a 
loss of coherence as momentum transfer to the nucleus increases. We use the Helm form factor 
parameterization F(E nr ) = 3ji(qr n )/qr n ■ exp(— (qs) 2 /2). We take the effective nuclear radius 

1.2A 1 / 3 fm and the skin thickness s 



T n = \Aq - 5s 2 , with r 
transfer to the nucleus is just q 



1 fm 
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The momentum 



-\/2MjvE„ r . 

We assume the dark matter halo velocity distribution to be Maxwellian, and perform the 
integration over f(v)/v following [30 . We take the velocity dispersion to be vq = 230 km s _1 , and 
as in 31 set the rotational speed of the local standard of rest v rot = vq. Considering the sun's 
velocity [BCJ implies an earth-halo velocity of 243 km s _1 . The lower limit of the integral in Eq. 6.1 
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Figure 9: For both XENON10 and XENON100 we show exclusion limits in the m x — a n plane, for each 
of the three C e ff curves shown in Fig. w (with corresponding line style). The shaded ellipse corresponds 
to a combined CoGeNT+DAMA 99% allowed region fct . Parameter space above the curves is excluded at 
90% C.L. 



is Vmin = \/l/(2mjvE nr )(m./vE m .//x). The upper limit is set by the galactic escape velocity, which 
we take to be v esc = 600 km s _1 |B2l. 



In calculating 90% C.L. exclusion limits, we follow the procedure described in Sec. 4.1.2. The 



combined low-energy efficiency r\ (Fig. |] and Fig. ^, 0) is applied directly to the predicted nuclear 
recoil energy spectrum (Eq. fO]). The spectrum is then convolved with the detector resolution as 
obtained directly from the simulation [Em. The stated cut acceptances are then applied. These 
are given for XENON10 by Table I of fu|, and XENON100 by Fig. 2 of @. For XENON100, 
we presume their total acceptance curve |l] already includes their estimate of ?7(S1). As shown in 
Fig. [| , this should roll off rapidly below 4 photoelectrons. Since we have already applied rj to 



Eq. 6.1, we take the stated XENON100 acceptance to be a constant 0.65 below 5 photoelectrons. 



We also note that the XENON100 acceptance box is defined between /i and the 15 electron S2 
software threshold (Tj. This may result in a slightly lower acceptance at low energies, relative to an 
acceptance box defined between /i and /i — 3<r. We assume such an effect, if present, is negligible. 



The lower limit of integration for calculating the number of events predicted by Eq. 6.1 



determined by the SI threshold, converted to keVr via £ e //- For both experiments, we use 2 pho- 
toelectrons for this lower limit. This is reasonable despite the stated 4 photoelectron analysis 
threshold of XENON100, as they observe no events above either threshold. If we had instead used 
a 4 photoelectron SI threshold for XENON100, those exclusion limit curves would have been about 
30% higher for m x < 10 GeV. We reiterate that £ e ff has only been extrapolated as shown in Fig. 
H, which was necessary to reproduce the observed XENON10 nuclear recoil band. We also note 
that because of the detector acceptance (77), any extrapolation below 2 keV is largely irrelevant. 

7. Summary 

Our analysis has shown that reasonable bounds on the energy dependence of C e ff can be deduced 
from existing data on Q v , in conjunction with nuclear recoil band data. These bounds, with their 
requisite extrapolation, are shown in Fig. |8|. This result makes use of the simple observation that 
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the energy dependence of N e /N 7 , and therefore also Q y /C e ff, are encoded in the nuclear recoil 
band logio(S2/Sl). It has already been pointed out that the sensitivity of liquid xenon detectors 
to light dark matter is substantially enhanced by considering the finite resolution Jy. However, we 
have shown that the detector resolution cannot be correctly described by only considering Poisson 
fluctuations in the number of detected SI photoelectrons. This is due in large part to the effect 



of taking 50% acceptance for nuclear recoils, as described in Sec. |L2j. We have made available 
pq] the discrete resolution distributions for each C e ff an d Q y case shown in Fig. |[ We have 
also shown a rigorous treatment of liquid xenon detector acceptance, for low energy nuclear recoils. 
The acceptance shown in Fig. [| and Fig. assume that position-dependence of signals leads to 
sub-dominant effects. If this were not true, the actual detector acceptance could be lower. 

For equal SI thresholds of 2 photoelectrons, XENON10 shows more sensitivity than does 
XENON100 to dark matter with particle mass m x < 30 GeV. This is a result of the difference 
in rj between the two experiments (compare Fig. with Fig. ffl), which in turn is driven almost 
exclusively by the higher L y /S e obtained by XENON10. An SI coincidence requirement of n > 3, 
which was not employed by cither experiment, would be a good example of another factor which 
that would significantly lower 77. Our results indicate that at 90% C.L. the XENON10 data are 
inconsistent with the combined CoGeNT+DAMA 99% C.L. allowed region p3|, given our central 
(Fig. 0, solid curve) assumption about £ e ff- With the most conservative assumption about C e ff 
(Fig. g, dashed curve), XENON10 (XENON100) is marginally (comfortably) consistent with the 
combined CoGeNT+DAMA 99% C.L. allowed region. Taken individually, either the CoGeNT or 
DAMA allowed regions |3^| remain consistent with the XENON10 constraints; these cases are not 
shown in Fig. |9| A separate analysis of the CoGeNT data in conjunction with a 50% exponential 
background contribution finds consistent parameter space for m x < 6 GeV |plj . We also note that 



our exclusion limits appear to be less stringent than those derived in 27 , for m x < 6 GeV. This is 
probably a result of the treatment of rj and the detector resolution. 

We note that neither astrophysical uncertainties nor uncertainty in the form factor F(E nr ) 
have been taken into account in the exclusion limits. Considering these factors will likely ease the 
tension between the null results of XENON10 and XENON100, and the dark matter interpretation 
of CoGeNT and DAMA. We also remind the reader that the consistency of Case 2 described in Sec. 



3.2 required an approximately la systematic decrease in the measured value of Q y from H. This is 
consistent with expectations, and it highlights the need to obtain both Q y and nuclear recoil band 
data under exactly the same operating conditions. We do not expect the acquisition of such data 
to significantly modify our conclusions, but rather to decrease the range of reasonable values of the 
scintillation yield from that shown in Fig. pi 
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A. Scintillation decay time constant r 

In this Appendix we briefly discuss the scintillation decay time constant (r) for liquid xenon. 
The relevance of this discussion centers on the choice of the XENON100 collaboration to require 
n > 2 coincidence in a 20 ns time window. For the 300 ns coincidence window used by XENON10 
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it is largely irrelevant, since > 99.9% of scintillation occurs within 300 ns in any case. It is well 
established that the excited molecular state Xe^ exists in either a triplet or singlet spin configuration, 
with a different relaxation time in each case |17j]. In other words, the primary scintillation signal 
as a function of time is 

(A.I) 



Sl(t) = A s • e (-*/^) + A t ■ e { - t/Tt) . 



Measurements of r s = 2 ns and Tt — 27 ns for electron recoils, and r s = 4 ns and Tt = 22 ns for 
nuclear recoils can be found in |Q. We are not aware of any physical reason why the singlet and 
triplet state decay times should depend on the incident particle species, and so assume t s = 4 ns 
and r t — 27 ns. 
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Figure 10: Idealized primary scintillation pulse shapes (solid curves) for nuclear (A s — 5.00) and electronic 
(A s — 0.35) recoils. In both case, A t = 1. The cumulative sum of each pulse is also indicated (dashed 
curves). Note that the pulse shapes, which are normalized to unit area, have been scaled by x 10 for display. 

The intensity ratio I s /It = A s T s /A t T t , is often used to indicate the relative strength of each 
component. For electron recoils, I s /It = 0.05 |l7j and for recoils resulting from alpha particles, 
I s /It = 0.43 [|34[ . We are not aware of any measurement for low-energy nuclear recoils from 
neutrons, but we can infer the intensity ratio from the prompt fraction of scintillation reported by 
the XENON10 collaboration (Fig. 2 of (ul). We calculate the prompt fraction f p as defined in 
p5| , H) for the idealized pulse defined by Eq. A.l, with no attempt to account for pulse shaping due 
to electronics. We choose A s = 0.35 and A t — 1 for electron recoils, since these values reproduce 
the measured I s /It = 0.05. This leads to f p — 0.19, in reasonable agreement with pq] . Noting 
that |lq] report an asymptotic value of f p for nuclear recoils that is about 12% larger than their 
electron recoil value, we then ask what value of A s this implies for nuclear recoils. It turns out that 
A s = 5.00 and A t — 1 (I s /It — 0.74) provide consistency with that work. 

This implies that about 75% of scintillation photons from nuclear recoils would be expected to 
appear in the first 20 ns, as summarized in Fig. nOl Note that if it were the case that A t = for 
nuclear recoils, 99% of scintillation photons would be expected in the first 20 ns. However, the data 
|17| , p4| do not support such a conjecture, nor does the prompt fraction reported in p5l B5| . Given 
these concerns, a direct measurement of the intensity ratio for nuclear recoils is probably in order. 
Whether or not that occurs, the XENON100 collaboration may wish to increase their coincidence 
time window to e.g. 100 ns, for a more comfortable 99% collection of the SI signal. 
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